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<N : Abstract 



For many years it was believed that an unstable periodic orbit with an odd number of real Floquet 



O 

■ multipliers greater than unity cannot be stabilized by the time-delayed feedback control mechanism 



of Pyragus. A recent paper by Fiedler et al. 



151 ] uses the normal form of a subcritical Hopf 



CD ' bifurcation to give a counterexample to this theorem. Using the Lorenz equations as an example, 

u' 

we demonstrate that the stabilization mechanism identified by Fiedler et al. for the Hopf normal 
form can also apply to unstable periodic orbits created by subcritical Hopf bifurcations in higher- 
dimensional dynamical systems. Our analysis focuses on a particular co dimension- two bifurcation 

(N 

' that captures the stabilization mechanism in the Hopf normal form example, and we show that 

l> 

the same codimension-two bifurcation is present in the Lorenz equations with appropriately chosen 
Pyragus-type time-delayed feedback. This example suggests a possible strategy for choosing the 
feedback gain matrix in Pyragus control of unstable periodic orbits that arise from a subcritical 
Hopf bifurcation of a stable equilibrium. In particular, our choice of feedback gain matrix is 



informed by the Fiedler et al. example, and it works over a broad range of parameters, despite 
the fact that a center-manifold reduction of the higher-dimensional problem does not lead to their 
model problem. 
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I. INTRODUCTION 



Time-delayed feedback control has been used as a method of stabilizing unstable periodic 
orbits (UPOs) or spatially extended patterns by a number of authors. The method of Pyra- 
gus [l|, sometimes called 'time-delayed autosynchronization' (TDAS), has attracted much 
attention. Here, the feedback F is proportional to the difference between the current and a 
past state of the system. That is, F = K(x(t — r) — x(t)) where x(t) is some state vector, 
t is the period of the targeted UPO and K is a feedback gain matrix. Advantages of this 
method include the following. First, since the feedback vanishes on any orbit with period 
r, the targeted UPO is still a solution of the system with feedback. Control is therefore 
achieved in a non-invasive manner. Second, the only information required a priori is the 
period r of the target UPO, rather than a detailed knowledge of the profile of the orbit, or 
even any knowledge of the form of the original ODEs, which may be useful in experimental 
setups. The method Jias been implemented successfully in a variety of laboratory experi- 

6], and chemical 3,181] systems, as well as in 



laser \S , plasma 



10. 



11 



121 ]; more examples can be found in a recent review by 



ments on electronic 
pattern-forming systems 
Pyragus 113 1. 

A paper of Nakajima [14J gave a supposed restriction on the method of Pyragus. It was 
believed that if a UPO in a system with no feedback had an odd number of real Floquet 
multipliers greater than unity, then there was no choice of the feedback gain matrix K for 



which the method o 
of Fiedler et al. 



Pyragus could be used to stabilize the UPO. However, a recent paper 
15l | gives a counterexample to this restriction. They add Pyragus-type 
feedback to the normal form of a subcritical Hopf bifurcation and show that the subcritical 
periodic orbit can be stabilized for some values of the feedback gain matrix. The Hopf 
normal form is two-dimensional, so the subcritical orbit has exactly one unstable Floquet 
multiplier. The mechanism for stabilizing the orbit is through a transcritical bifurcation 
with a stable delay-induced periodic orbit. Just et al. [16] investigate a series of bifurcations 
in this system, which has the attractive feature that, despite the presence of the delay terms, 
much of the analysis can be carried out analytically. 

The subcritical Hopf bifurcation of a stable equilibrium is a generic mechanism for cre- 
ating UPOs with an odd number of unstable Floquet multipliers. Such bifurcations occur 
in a number of physical systems, such as the Belousov-Zhabotinsky reaction-diffusion equa- 



2 



tion 
lasers 



17L the Hodgkin-Huxley model of action potentials in neurons 18j, and in NMR 
191 ] . The reduction of these higher- dimensional dynamical systems to the two- 
dimensional normal form of the Hopf bifurcation problem is a standard procedure 2C|, 121 ]. 
Moreover, if Pyragus-type feedback delay terms were added to the model ODEs, then these 
(infinite-dimensional) dynamical systems could likewise be reduced to the standard two- 
dimensional normal form in a vicinity of a Hopf bifurcation 22|, with the parameters of 
the feedback control matrix K modifying the coefficients in the normal form. Despite this 
disconnect between center manifold reduction of delay equations to Hopf normal form, and 
the example of Fiedler et al. in which the feedback delay terms are added directly to the 
Hopf normal form, we find that the same stabilization mechanism of subcritical Hopf orbits 
applies to both their example and to the one we present for the Lorenz equations. 

Specifically, we study a subcritical Hopf bifurcation of a stable equilibrium in the Lorenz 
equations [23l. |24| . and show that Pyragus-type feedback can stabilize the subcritical periodic 



orbit. As in the example in 1^], in the absence of feedback, the bifurcating periodic orbit 
has exactly one real unstable Floquet multiplier. It also has one stable Floquet multiplier, 
and one Floquet multiplier equal to one (corresponding to the neutral direction along the 
orbit). The 3x3 gain matrix multiplying the Pyragus feedback terms can be chosen in 
many different ways. We give two examples in which we choose the structure of the gain 
matrix in different ways and show that they give quite different results. 

n our first example, we choose the gain matrix in a manner suggested by the results 



m 



3 



and 



161 ]: there is no feedback in the stable direction of the UPO, and Pyragus- 



type feedback in the direction of the unstable Floquet multiplier, which is identical in form 
to the feedback in [15j. In this way, the problem of choosing the nine parameters in the 
3x3 gain matrix is reduced to one of making an informed choice of the two parameters 
employed in [2]. We find that the subcritical orbit can be stabilized over a wide range 
of values of our two bifurcation parameters: the amplitude of the feedback gain, and the 
usual control parameter p in the Lorenz equations. We identify a codimension-two point 
in the Hopf normal form example, where two Hopf bifurcations collide, and show that 
the same codimension-two point can be found in the Lorenz system with this choice of 
feedback, and the bifurcation structure is qualitatively the same in the two cases. This 
codimension-two point captures the stabilization mechanism in both examples: the periodic 
orbits created by the two Hopf bifurcations exchange stability in a transcritical bifurcation. 
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The curve of transcritical bifurcations in our two-parameter plane emanates from the Hopf- 
Hopf codimension-two point. 

Our second choice of the gain matrix is a real multiple of the identity. This is also a 
natural choice, but in contrast with our first example we show that here the subcritical 
orbit cannot be stabilized for any parameters close to the original Hopf bifurcation. Our 
two bifurcation parameters are again the amplitude of the gain and the parameter p in the 
Lorenz equations. We give analytical results on the location of Hopf bifurcation curves and 
hence deduce the stability of the periodic orbit as it bifurcates. 

This paper is organized as follows. In section [III we review some results from Fiedler et al. 
[15I and Just et al. [lfj . In section II III we give our example system of the Lorenz equations 
with Pyragus feedback. We give two examples of the choice of gain matrix. We explain for 
the first example how we choose the gain matrix to stabilize the subcritical Hopf orbit, and 
show that the bifurcation structure of this system is the same as that for the normal form 
system. For the second example, the gain matrix is a real multiple of the identity and we 
show that the subcritical orbit cannot be stabilized. Section HVl concludes. 

II. THE HOPF NORMAL FORM WITH DELAY 

In this section we recap the results of [ijj] and identify a particular codimension-two point 
in the Hopf normal form with delay which we will later examine for the Lorenz equations 
with feedback. This codimension-two point acts as an organizing center for the bifurcations 
involved in the mechanism for stabilizing the periodic orbit. 

The normal form of a subcritical Hopf bifurcation with a Pyragus-type delay term is: 



with z G C, and parameters A, 7 G R. The feedback gain b = b^e 113 G C, and the delay r > 0. 
The linear Hopf frequency has been normalized to unity by an appropriate scaling of time. 
We consider A as the primary bifurcation parameter. We consider only 7 < 0; this is the 
case in the Lorenz example. 

For the system with no feedback (i.e. b = 0) we can write z = re w and then 



z{t) = (A + i)z{t) + (1 + i-f)\z(t)\ 2 z(t) + b{z{t - t) - z{t)) 



(1) 




(2) 
(3) 



= 1 + 7 r 2 . 
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Periodic orbits exist with amplitude r 2 = —A if A < 0, so 9 = 1 — 7A and the orbits have 
minimal period T = 2tt/(1 — 7A). We refer to these orbits as the Pyragus orbits, and it is 
these orbits that we wish to stabilize non-invasively by adding an appropriate feedback term 
(i.e. with b ^ 0). 



Following 



15l ] and [if]], we define the Pyragus curve t = Tp(X) in A-r space, along which 



the feedback vanishes on the Pyragus orbits: 

^ = T^V M 

We plot this curve in A-r space in figure dj along with curves of Hopf bifurcations from the 
zero solution. In later sections, we set r = Tp(A), as our main purpose is the non-invasive 
stabilization of the Pyragus orbits. 

The zero solution of ([I]) undergoes Hopf bifurcations when the characteristic equation 
has purely imaginary solutions. Setting z(t) = e vt in ([1]) and linearizing we find: 

v = X + i + b(e rrr - 1). 

Writing rj = iu and separating into real and imaginary parts gives 

= A + & [cos(/3 - urr) - cos/3], (5) 
uj — 1 = & [sin(/3 — ojt) — sin/3]. (6) 

These equations define the Hopf curves r = r#(A), in A-r space, parameterized by the linear 
frequency uj associated with the bifurcating periodic orbit. There are multiple branches to 
this curve, which we show in figure [0(a), but we concentrate on the one which intersects the 
curve r = rp(A) at (A, r) = (0, 2tt). The solution of the characteristic equation at A = 0, 
r = 2n has u = 1 and corresponds to the Hopf bifurcation to the Pyragus orbit. 

Figure [H shows the possible configurations of the curves r = rp(A) and r = r#(A) as the 
parameter b is varied. The curves typically cross in two places: at A = 0, and at a second 
location depending on 6q- 

At b = 6q, the two curves are tangent at A = and only intersect 



once. Just et al. 



161 ] show that 



27r(7 sin (3 + cos 0) 
For simplicity, we assume 6g > 0, so we must have 7sin/3 + cos/5 < 0. 
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(a)6 = 0.1 (b)6 = 0.025 (c)6 = 0.0214 

FIG. 1: The figures show the curves r = Tp(A) (dashed curve) and r = Tff(A) (solid curve) for 
three values of 6q. The remaining parameters in ([TJ) are (3 = tt/4 and 7 = —10. In (a), portions of 
four of the Hopf curves are shown, but in (b) and (c) we show only the curve which passes through 
(A,r) = (0, 27r). In all three cases, the two curves cross at A = (shown by a solid dot). In (a), 
(with bo > &q), the curves cross again in A < (shown by an empty dot), and in (c), (60 < 6q) the 
curves cross again in A > 0. Case (b) has bo = 6q and the two curves are tangent at A = 0. The 
origin is stable (unstable) in those regions marked by an s (u). 

We define a curve of Hopf bifurcations b = &Q° pf (A) in X-bo space by the location of 
the second intersection of Tp(A) and Th(A). This is a Hopf bifurcation to a delay-induced 
periodic orbit, that is, a periodic orbit arising from the addition of the delay terms; one for 
which the feedback does not vanish. 

We adopt the convention that a Hopf bifurcation of a stable equilibrium is called 'super- 
critical' ('subcritical') if the resulting periodic orbit bifurcates into the parameter regime 
where it coexists with the unstable (stable) equilibrium. Such a supercritical bifurcation 
generically produces a stable periodic orbit 20], while the subcritical case produces an un- 
stable periodic orbit. In the absence of feedback, the bifurcating orbit is subcritical and 
unstable. The mechanism for stabilization involves the additional delay-induced Hopf bi- 
furcation at 6^ opf (A). This bifurcation can change the trivial equilibrium from being stable 
to being unstable. Consequently, the Pyragus orbit may then co-exist with an unstable 
periodic orbit. As the Hopf bifurcation at 6^ opf (A) passes through A = 0, the original Hopf 
bifurcation to the Pyragus orbit at A = changes from a subcritical one to a supercritical 
one. This is all done without otherwise altering the form of the Pyragus orbit. 

The Hopf bifurcation of the zero solution to the Pyragus orbit at (A, r) = (0, 2tt) changes 
from subcritical to supercritical as described above as bo is increased through b$, since for 
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bo > 6g, the curve rp(A) lies 'inside' r#(A). In this sense, 6q is the smallest value of the 
feedback gain for which the Pyragus orbit is stabilized immediately after the bifurcation 
point. The minimum positive 6q can be selected by choosing ft such that 7 = tan ft. 



A. A codimension-two bifurcation point 

We now review some of the details of the bifurcation structure of the system (CQ) which 
are described in Just et al. [la ], and identify the codimension-two point we examine in the 
Lorenz system. We consider A and b as two bifurcation parameters, and fix r = Tp(A). 

The mechanism by which the Pyragus orbit is stabilized is through a transcritical bifur- 
cation with a delay-induced periodic orbit. As shown in Just et al. the transcritical 
bifurcations occur when 

-1 

T — 

6 (cos ft + 7sin/3) ' 
or, in A-&o space, since r = rp(A), when 

A = -(1 + 2tt& (cos/3 + 7 sin/?)) = - (l - ^ 
7 7 V b c 



This line of transcritical bifurcations collides in A-60 space with the two curves of Hopf 
bifurcations A = and b = &Q° pf (A) at (A, 60) = (0, b^), at a double-Hopf codimension- 
two point. In figure [2] we sketch the bifurcation structure around this point in A-60 space. 
From this figure we can see that in order for the Pyragus orbit to bifurcate stably (i.e. 
supercritically) at A = 0, we must have 60 > b^. 



III. THE LORENZ EQUATIONS WITH TIME-DELAYED FEEDBACK 

We now use the Lorenz equations as an example system to demonstrate that the feedback 
described above can also stabilize orbits arising in a subcritical Hopf bifurcation in a higher- 
dimensional system of differential equations. We give two examples of a choice of feedback 
gain matrix. The first choice is informed by the results given above, and for the second 
choice we set the gain matrix equal to a real multiple of the identity. In the first example, 
we further locate the codimension-two point described in section III in the Lorenz system 
with feedback, and show that the bifurcation structure is the same as in the normal form 
case. 
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FIG. 2: In (a), the two Hopf bifurcation curves A = and bo = & Hopf (A) (solid bold lines), and the 
transcritical (TC) bifurcation curve (dashed bold line) divide the X-bo plane into five regions. The 
curves intersect at (A,6q) = (0, 6q). The Pyragus orbit is stable in the shaded region, (b) shows a 
schematic representation of the solutions as a path C is traversed anticlockwise around the origin. 
Solid lines represent stable solutions and dashed lines represent unstable solutions. 



The Lorenz equations 23j, |24j are most often written in the following form: 



x = a(y - x), 
y = px — y — xz, 
z = —az + xy, 

for real parameters a, a and p. Lorenz and most other authors studied the parameter 
regime a = 10, a = 8/3, p > 0, and we continue in the same manner. Taking p as the 
primary bifurcation parameter, the zero solution is stable for p < 1 and loses stability in a 
supercritical pitchfork bifurcation at p — 1. Two further equilibria are created at 



x± 



(±VMp-i)), i^Wp-i)), p - 1). 



As p is increased further, these equilibria each undergo a subcritical Hopf bifurcation at 



Ph 



a (a + a + 3)/(cr - a - 1) w 24.74 (see |20| and [2l| for further details). It is this 
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FIG. 3: The top figure is a bifurcation diagram of the subcritical Hopf bifurcation in (|T|), showing 
the fixed point at zero and the bifurcating branch of unstable periodic orbits. Solid lines indicate 
stable solutions and dashed lines indicate unstable solutions. The lower figure shows the period, 
T, of the bifurcating periodic orbit, as a function of p. 

bifurcation that we study in the following, so we shift coordinates to be centered around x + 
and rescale to obtain: 

u = a(v — u), 

v = u — v — (p — l)w — (p — l)uw, (7) 
w = a{u + v — w + uv). 

Figure [3] shows a bifurcation diagram of the subcritical bifurcation of the zero solution of (]7j), 
and also the period T of the bifurcating orbits, which we use to determine the delay time 
Tp{p) in the controlled system. The periodic orbit exists for 13.926 < p < ph] at the lower 
boundary it collides with a fixed point in a homoclinic bifurcation. 

A. Adding time-delayed feedback 

We now add Pyragus-type feedback to the Lorenz equations. We write 











' ' u T — 




V 


= J(p) 


V 


+ N(u,v,w) + r 


v T — V 


(8) 
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where 



/ 



J(P) 



—a a 
1 -1 

a a 





<P-l 
—a 



\ 



( 



I 







V 



(9) 



/ 



u T = u(t — r), etc. and T is a_3 x 3 rea 
section we use the results of 



N(u,v,w)= -(p-\) uw 
auv 

feedback gain matrix, to be determined. In this 
151 ] and 16| to inform our choice of the control matrix T. In 
general, T would contain nine independent parameters, but the method we describe reduces 
this to only two. In section [III CI we describe the dynamics when T is a real multiple of the 
identity. Note that if this system were reduced to normal form around the Hopf bifurcation 
point, the resulting equations would not be the same as ([I]). That is, there would be no 
delay terms; the delay terms here would only have the affect of altering the parameters in 
the usual Hopf normal form. See 22j for more details. 

At the bifurcation point (p = ph), J has one real negative eigenvalue (which we denote 
by —\h = —a — a — 1 m —11.7), and a pair of purely imaginary eigenvalues (iza^, Uh = 



\/ (2acr(a + l))/(er — a — 1) ~ 9.62). The center manifold of the original problem with no 
feedback is therefore two-dimensional, and the eigenvectors of J can be found explicitly 
(see {21]). Close to the bifurcation point, the subcritical orbit will lie in a two-dimensional 
manifold which is close to the center subspace at the bifurcation point. We therefore choose 

r = ege-\ 

where E is the matrix of eigenvectors which puts J(ph) = Jh i n Jordan normal form, that is 



( -\ h ^ 



E~ l J h E 



\ 



-u h 
u h 



(10) 



and 



/r 



no 

G= 6 cos (3 —b sin j3 \ (11) 
\0 bo sin p b cos (3 , 

There is then no feedback in the stable direction, and Pyragus-type feedback in the directions 
tangent to the center manifold. 

In the following numerical results, we set (3 = 7r/4, as in Fiedler et al, and vary bo- 
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(a)6 = 1.2 




(b)6o = 0.1 



FIG. 4: The figure shows the curves r = Tjj(p) (solid lines) and r = rp(p) (dashed lines) for 
&o = 1-2 > 6§ and 60 = 0.1 < 6q. Remaining parameters are a = 10, a = 8/3 and (5 = 7r/4. 
Compare with figure [TJ The Hopf bifurcation at p = p^ ~ 24.74 is shown with a solid dot. In (a) 
the curves additionally cross in p < p^ (shown by an empty dot). The origin is stable (unstable) 



in those regions marked with an s (u). These figures were produced using DDE-BIFTOOL 
B. Numerical results 



We use the continuation package DDE-BIFTOOL [25[ to analyze the delay-differential 
equation (jS]). The primary bifurcation parameter is p, with the subcritical Hopf bifurcation 
for the system without feedback occurring at p = p^ « 24.74. Recall we have set a = 10, 
a = 8/3 and (3 = tt/4. 

First, we locate Hopf bifurcations of the trivial solution in the p-r plane, for various 
values of bo. Figure H] shows curves of Hopf bifurcations r = Th(p) for 6 = 1-2 and b = 0.1. 
We also plot the curve r = r P (p), given by the period of the bifurcating subcritical orbits 
(see figure [3]) • Figure H] is qualitatively similar to figured] (the corresponding figure for the 
normal form case). For 60 = 1.2, Tp(p) lies inside th(p), and so we expect, by analogy with 
the normal form case, that choosing r = Tp(p) will stabilize the Pyragus orbit. For bo = 0.1, 
Tp(p) lies outside Th(p), and so the feedback cannot stabilize the Pyragus orbit near onset, 
since it bifurcates subcritically (i.e. it coexists with the stable equilibrium from which it 
bifurcates). The codimension-two point occurs at some value of bo that is the boundary 
between these cases. 

We use DDE-BIFTOOL to locate this codimension-two point. As in the normal form case, 
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FIG. 5: (a) shows curves of Hopf bifurcations from the zero solution in the p-bo plane for the 
Lorenz system with delay. The vertical line is the bifurcation to the Pyragus periodic orbits, and 
the second curve is the bifurcation to the delay- induced periodic orbits. The curve of transcritical 
bifurcations of periodic orbits is not shown, but it can be seen in figure El The dotted ellipse is 
the curve traversed to generate the bifurcation diagram in (b). Here, solid lines indicate stable 
solutions, and dashed lines indicate unstable solutions. Parameter values are a = 10, a = 8/3, 
P = 7r/4. The ellipse is parameterised by 9: p — ph = 1.4 cos 9, bo — b$ = 0.03 sin 9, and encloses the 
codimension-two point. This figure was generated using DDE-BIFTOOL. Compare with figure [2j 

we set t = Tp(p). We do not have an analytic form for Tp(p), so we numerically estimate 
Tp(p) in the following way. For p < ph we set rp(p) equal to the period of the bifurcating 
period orbits for the system with no feedback (see figure [3]). We want to continue Tp(p) 
into p > ph, so we can complete both sides of the bifurcation diagram, so here we set 
Tp (p) = Th /(1 - B(p - p h )), where B = -0.0528, and r h = r P {p h ) « 0.6528. This choice of 
B ensures that Tp(p) is continuous and has continuous first derivative at p = ph- 

With the parameter restriction r = Tp(p), we generate curves of Hopf bifurcations from 
the zero solution in the p-bo plane; these are shown in figure [5j We can then estimate the 
location of the codimension-two point, at (p, bo) = (ph, &o)> the point where the two curves of 
Hopf bifurcations cross. We find 6q ~ 0.221. We follow a path around the codimension-two 
point and track the amplitude and stability of the bifurcating periodic orbits; a bifurcation 
diagram of the periodic orbits is shown in figure [H The transcritical bifurcation of periodic 
orbits can clearly be seen. Note that figure [5] is qualitatively similar to figure El showing 
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FIG. 6: The figure shows a contour plot of the modulus of the largest Floquet multipliers for the 
Pyragus orbit as p and 60 are varied. The shaded area indicates the region where the periodic orbit 
is stable. The codimension-two point is marked with a dot. The boundary of the stable region 
which emanates from this point is the transcritical bifurcation of periodic orbits. The left boundary 
of the stable region corresponds to a period-doubling bifurcation. This figure was produced using 
DDE-BIFTOOL. 

that the bifurcation structure in the normal form case, and in our Lorenz example are the 
same. 

With the additional feedback, the Pyragus orbits are stable for a wide parameter range. 
In figure El we show the stability of the Pyragus orbits as p and 60 are varied. The tran- 
scritical bifurcation can be seen as the boundary of the stable region which terminates at 
the codimension-two point. The orbits also undergo an instability at around p ~ 17. Along 
this boundary of the stable region, the periodic orbits have a Floquet multiplier equal to 
— 1, and so the instability is a period-doubling bifurcation. 

In figure [7] we show results of forward time integration of the delay-differential equa- 
tion (jSJ), at p = 23 < ph, with r = r P (p) = 0.7191. Initally, b = 1.2, and the Pyragus orbit 
is stable. The feedback is then turned off (i.e. 60 = 0) at t — 50, and the trajectory decays 
back to the zero solution. We have used DDE-BIFTOOL to confirm the stability of these 
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FIG. 7: (Color online) The figure shows time integration of the Lorenz equations with feedback 
which is switched off at t = 50 (indicated by a dashed vertical line). Parameter values are a = 10, 
a = 8/3, (3 = tt/4, p = 23 and r = r P (p) = 0.7191. b = 1.2 for t < 50 and b = for t > 50. The 
Pyragus orbit is initially stable, but without feedback, the trajectory decays back to the origin. 
The data was computed using the Matlab routine dde23 for integrating delay-differential equations. 

orbits. 

The structure around the codimension-two point also tells us that the delay-induced orbits 
can be stable in the region p > p h . For example, at p = 24.8388, b = 0.22, r = 0.6494, we 
can use DDE-BIFTOOL to show that there exists a stable delay-induced periodic orbit with 
a period of 0.6537. Figure [8] shows time integration at these parameter values, with feedback 
turned on at t = 5. The figure shows the chaotic attractor for t < 5 and an approach to a 
stable periodic orbit for t > 5. 



C. A second example 

Another natural choice for the gain matrix T is a real multiple of the identity. As a 
comparison to the results given above, we now consider this case, so write 



/in n\ 



b 



1 
1 
1 



bo e 



(12) 



For this choice of T we can obtain analytical results, and we show that the bifurcation 
structure is different from our previous example, whatever the value of bo. We show first that 
unstable periodic orbits with real positive Floquet exponents cannot be stabilised using this 
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FIG. 8: (Colour online) Figures (a) and (b) show two time integrations of the Lorenz equations. 
Figure (a) has no feedback terms, and in figure (b) the feedback is 'turned on' at t = 5. Without 
feedback, the well-known chaotic strange attractor is stable, but with the feedback, a delay-induced 
periodic orbit becomes the stable solution. The stability of the periodic orbit has been confirmed 
using DDE-BIFTOOL. Parameters are p = 24.8388, b = 0.22, r = 0.6494. The stable periodic 
orbit has a period of 0.6537. Figure (c) shows the time plot of figure (b). The data was computed 
using the Matlab routine dde23. 

form of feedback. We then discuss the shape and location of the curves of Hopf bifurcation 
from the origin, in a similar manner to the previous section, to give a comparison of the 
two types of feedback. The co dimension-two point described previously does not exist, and 
the Hopf bifurcation to the Pyragus orbit is always subcritical. That is, the periodic orbit 
in p < ph always bifurcates unstably from an equilibrium which is stable in p < ph and 
unstable in p > p h . 



15 



Consider a Pyragus solution u*(t) of (JSJ), which is periodic with period r, and with V 
as in ffl2l) . Then if p is a Floquet exponent of u* (in the system with no feedback), the 
characteristic equation for u* in the system with feedback is given by 

(A-6 (e- Ar -l))=/i- (13) 

Both here, and in the analysis which follows below, simplification is possible because T is a 
multiple of the identity. If u* has one Floquet exponent p which is real and positive, then 
it can be shown (see e.g. 26J]) that there always exists at least one solution A which is real 
and positive. Hence stabilization of u* cannot be achieved. 

We note that at the subcritical Hopf bifurcation from the zero solution in ([7]), the orbit 
which bifurcates has one real stable multiplier (inherited from the zero solution) and one 
neutral multiplier, and so the remaining unstable multiplier must be real. Therefore the 
Pyragus orbit cannot be stabilized close to the Hopf bifurcation using this type of feedback. 

In addition, we now follow the method of the previous section to find curves of Hopf 
bifurcations of the zero solution of ([8]) in p-r space. This provides us with a comparison 
of the two types of feedback used in this and the previous section. We are particularly 
interested in those curves which pass through (p, r) = (p^r^), with Hopf frequency equal 
to u>h, as this is the Hopf bifurcation to the Pyragus orbit. 

Consider the linearisation of (JSj) about the origin, and write (u,v,w) T = ue xt . Then 

\ue xt = J{p)ue xt + Tu{e- XT - l)e xt , 

and so for a non-trivial solution (u ^ 0) to exist, A must satisfy the characteristic equation: 

det[(A - 6 (e- Ar - 1))/ - J(p)\ = 0. (14) 

Note that we have been able to simplify this equation because in this example T is a multiple 
of the identity. This tells us that g(X) = (A — bo(e~ Xr — 1)) are the eigenvalues of J(p). 

When p is close to ph, J(p) has one negative eigenvalue — Ai(p), and a complex conjugate 
pair we denote as p(p) ± iv(p). Note that p(p^) = 0, and v(p h ) = u h . 

We find curves of Hopf bifurcations in p-r space, by writing A = iuo (u e IR) and setting 
giiui) equal to the eigenvalues of J(p), p(p) + iv(p). We note that setting g{iu) equal to the 
negative eigenvalue of J(p), — Ai(p), does not produce any Hopf curves which pass through 
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FIG. 9: The figure shows curves of Hopf bifurcations from the zero solution of ([8]) with V = bol. 
Curves are shown for bo = 0.1,0.05 and —0.1 (from right to left). The zero solution is stable to 
the left of these curves and unstable to the right of these curves. The dashed curve is the curve 
r = Tp{p). For each value of bo, the curve lies between p = p^ and p = p*(bo). 

(p, r) = (ph, Th), so we do not consider these here. Equating real and imaginary parts gives 

-6 (cos(^r) - 1) = p(p), (15) 
uj + bo sin(cur) = u(p). (16) 

Equations (fl5j) and ffTB"]) describe curves of Hopf bifurcations in p-r space, parameterized by 
the Hopf frequency uj. Figure [9] shows examples of the shape of these Hopf curves; compare 
with figure H] showing the Hopf curves in our previous example. The zero solution is unstable 
to the right of the curves and stable to the left of the curves. We now explain why we expect 
the curves to have this shape. 
We have: 

(v(p)-uj) 2 = p(p)(2bo-p(p)). 

Note that we need /z(p)(2& — £*(/?)) > for solutions to exist. That is, if bo > 0, we need 
< n(p) < 26 , and if bo < 0, we need 26 < p(p) < 0. Since p(p) is a monotonically 
increasing function over the range of p we are considering, this gives a connected range of 
p for which Hopf bifurcations can occur, with boundaries at p = ph (since p{pu) = 0) and 
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at p = p*(&o) (where p(p*) = 2b Q ). Note that at p = ph, oj = v(ph) = Uh, and so this is the 
Hopf bifurcation to the Pyragus orbit. 

The solutions for r along the Hopf curve solve 



cos(c<;t) = 1 — 




Since p(ph) = 0, at p = ph, ojt = 2irn, and similarly since p(p*) = 2b , at p = p*, 
lot = 7i(2n + 1). The curve of Hopf bifurcations is thus tangent to the lines p = ph and 
p = p*, and forms a series of wiggles between these value of p. In particular, the curve is 
tangent to p = ph at (p, r) = (p h ,r h ), where u = u h , the Hopf bifurcation to the Pyragus 
orbit. 

Therefore, for any value of bo, the Pyragus curve r = rp(p) (which originates at (p, r) = 
(p/ij T h), and is also shown in figure E]), will always be to the left of the Hopf curves at p = ph- 
The zero equilibrium is therefore stable along the Pyragus curve close to the bifurcation point 
in p < ph-, and so the Pyragus orbit will always bifurcate unstably. Figure fTUl shows a contour 
plot of the largest Floquet multipliers of the periodic orbit, as p and bo are varied. For all 
values shown the periodic orbit is unstable. 



IV. DISCUSSION 

In this paper, we have demonstrated how the mechanism used by Fiedler et al. [ijj] 
for stabilizing periodic orbits with one unstable Floquet multiplier carries over to higher 
dimensional systems, using the Lorenz equations as an example. We use the results from 
their idealized example to inform our choice of the feedback gain matrix, and this method 
follows a set prescription which we expect could be used on other systems. First we find 
the two-dimensional linear center eigenspace of the system with no feedback at the Hopf 
bifurcation point. Feedback is then added in the directions lying tangent to this center 
subspace, using a 2 x 2 gain matrix of the form given by Fiedler et al. This then leaves 
only the two parameters (3 and bo to be chosen. For the example of the Lorenz equation, 
the subcritical orbits are stabilized over a wide range of parameters. We contrast this with 
an example of choosing the gain matrix as a real multiple of the identity. In this case the 
Pyragus orbit could not be stabilised. 

Choosing the gain matrix for our Lorenz equations example required a knowledge of the 
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FIG. 10: The figure shows a contour plot of the magnitude of the largest Floquet multipliers for 
the Pyragus periodic orbit, as p and 60 are varied. The periodic orbit is unstable for all parameters 
shown. This figure was produced using DDE-BIFTOOL. 

linearization of the system at the bifurcation point. This method may also be applicable 
in systems for which the governing equations are not known, if it is possible to access 
perturbations to the equilibrium solution near the Hopf bifurcation point, and hence extract 
the unstable eigenvectors numerically from experimental data. 

We have additionally shown that the Lorenz equations example contains a codimension- 



two point which is also present in the normal form example of 15|. This double-Hopf point 
is not generic. In the normal form example of [ljj] there is an additional SO (2) symmetry 
which is not present in the Lorenz example. Additional structures in the problem force a 



normally codimension-three phenomena 27| to be codimension-two, since the frequencies of 
the bifurcating periodic orbits are forced to be in one:one resonance at the codimension-two 
point. It would be of interest to examine this degeneracy in more detail, by understanding 
the mathematics behind the structure of the Hopf- Hopf bifurcation in these examples. We 
intend to investigate further examples to see how robust this bifurcation structure is, for 
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example, whether it appears in say, the Hodgkin-Huxley [18)] or Belousov-Zhabotinsky 17 ] 
examples. 
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